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SUMMARY 

In this paper, we establish the global well-posedness of the Cauchy problem for the Gross- 
Pitaevskii equation with an rotational angular momentum term in the space M^. 
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1 Introduction 

The Gross-Pitaevskii equation (GPE), derived independently by Gross [8] and Pitaevskii [15], 
arises in various models of nonlinear physical phenomena. This is a Schrodinger-type equation 
with an external field potential Vextit,x) and a local cubic nonlinearity: 

ihdtu+ —Au = Vextu + l3\u\^ u. (1.1) 
2m 

The GPE (jl.ip in physical dimensions (2 and 3 dimensions) is used in the meanfield quan- 
tum theory of Bose-Einstein condensate (BEG) formed by ultracold bosonic coherent atomic 
ensembles. Recently, several research groups [9, 12-14] have produced quantized vortices in 
trapped BECs, and a typical method they used is to impose a laser beam on the magnetic 
trap to create a harmonic anisotropic rotating trapping potential. The properties of BEC in a 
rotational frame at temperature T being much smaller than the critical condensation temper- 
ature Tc [10] are well described by the macroscopic wave function u(t,x), whose evolution is 
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governed by a self-consistent, mean field nonlinear Schrodinger equation (NLS) in a rotational 
frame, also known as the Gross-Pitaevskii equation with an angular momentum rotation term: 

^2 

ihdtu + — Am = V(x)u + NUq \u\^ u - nL^u, x G M^, t > 0, (1.2) 
2m 

where the wave function u{t, x) corresponds to a condensate state, m is the atomic mass, H is 
the Planck constant, N is the number of atoms in the condensate, O is the angular velocity 
of the rotating laser beam, and V{x) is an external trapping potential. When a harmonic 
trap potential is considered, V{x) — ^ {^i^i +^2^2 +^3^3) '^ith ui, L02 and being the 
trap frequencies in the xi-, X2- and xa-direction, respectively. The local nonlinearity term 
NUo\uf' u arises from an assumption about the delta-shape interatomic potential. Uq = 
Ai:h^as/m describes the interaction between atoms in the condensate with (positive for 
repulsive interaction and negative for attractive interaction) the s-wave scattering length, and 
Lz — —ih{xidx2 — X2dxi) is the third component of the angular momentum L = x x P with 
the momentum operator P = —ihV. 

After normalization, proper nondimensionalization and dimension reduction in certain 
limiting trapping frequency regime [16], it turns to be the dimensionless GPE in d-dimensions 
(d=2,3): 

1 9 ^ 

iut + -Au = Vd{x)u + j3d \u\ u - ni^u, x eM.'^,t>0, (1.3) 

where = i(xidx2 — X2dxi) and 



/3yW2^, ilWi+llA)/'^. d = 2, 

(3, ^ I hWi+jixl+jlxl)/2, d = 3, 



with 71 > 0, 72 > and 73 > constants, /3 = ao = y and = min{wi, W2, ^3}- 

In general, it is a rather complicated process about the dynamics of solutions (in particu- 
lar, vortex) for GPE (|1.2p under the interaction of trapping frequencies and angular rotating 
motion. The recent numerical simulation of GPE ()1.2|1 for different choice of trap frequencies 
(71, 72) can help us to understand the complicated dynamical phenomena caused by the angu- 
lar rotating and spatial high frequency motion. The case of different frequency 71 7^ 72 gives 
much complicated behavior and thus is rather difficult to be studied rigorously [2,3]. To our 
knowledge, the equation (|1.2p has been only investigated for some specific cases by numerical 
simulation. Therefore, to develop methods for constructing analytical solutions of the GPE 
l|l.l|) or some specific cases is the first step in order to understand the dynamics caused by the 
trapping and rotation. 

To begin with, we first consider the case 71 = 72 = w which means the spatial isotropic 
motion and focus on the Cauchy problem of the Gross-Pitaevskii equation with an angular 
momentum rotational term in two dimensions 

iut + -Au = Y \x\'^ u + (3 \u\^'^ u - LoL^u, a; e K^ t > 0, (1.5) 
u(0,x) = ito(a;), X e (1.6) 

where the wave function u = u{t,x) : [0, 00) x ^ C corresponds to a condensate state, 
A is the Laplace operator on R^, uj > 0, jS > and a G [l/u;,oo) are constants, and = 
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—i{xidx2 — X2dxi) = i{x2dxi — a^i^^a) is the dimensionless angular momentum rotational 
term [2,3]. We assume that the initial value 

uo{x) e^:^{ue H^R^) : \x\ue L^(M^)}, (1.7) 

with the norm 

It is clear that p.Sp is a special case of (|1.3p with a spatial isotropic trapping frequency 
(i.e. 7i =72) and f2 = a; in 2-dimensions. Note that for multi-dimensional GPE (|1.2|) . nothing 
is known about the exact integration except for the case when 71 = 72 73 for 3D) without 
the angular momentum rotational term, namely, $1 = considered in [5,6]. 

There are three ingredients that play important roles in the proof of our result. The first 
involves the solution of the Cauchy problem to the linear equation 

1 

iut -\- -Am = — |x| u - ujL^u, a; G M^, t ^ 0, 

Z Z (,l-"j 

u{0,x) ~ uo{x), x G M^, 

which is significant for investigating the properties of the evolution operator corresponding 
to the linear operator z^t + iA — ^|a;|^-|- uL^- The second one is to obtain the Strichartz 
estimates for the foregoing linear operator. The last one is that there exist two Galilean 
operators J{t) and H{t) (as blow) which can commute with the linear operator and can be 
viewed as the substitute of V and x respectively in the non-potential case. 
Now we state the main result of this paper. 

Theorem 1.1. Let uq G ^ o,nd p G [2, 00). Then, there exists a unique solution u{t,x) to the 
Cauchy problem (jl.5p - (|1.6p . And the solution satisfies, for any T G (0, 00), that 

u{t,x),J{t)u{t,x),H{t)u(t,x) e C{R; L^iR^)) n r'^-P^ {0,T; LfiR^)), V< G [0,T], 

where ~ 5 ^ p- '^i^) '^"■'^ -^(^) '^'^s defined as below as in ()2.7p and (|2.8p . respectively. 

Remark 1.2. Since the GPE ()1.5p (or ()1.3p ) in a rotational frame is time reversible and 
time transverse invariant, the above result is also valid for the case when t <0. 

The paper is organized as follows. In Sec. [21 the evolution operator of the linear equation 
and the Strichartz estimates about the former operator are first established. Sec. [3] is devoted 
to the derivation of some conservation identities such as the mass, the energy, the angular 
momentum expectation, the pseudo-conformal conservation laws in the whole space for 
(|1.5P " (|1.6p . Finally, the nonlinear estimates and the proof of Theorem 11.11 are obtained in 
SecEl 

2 The Strichartz estimates and some main operators 

Let u{t) be the solution of the linear equation (jl.Sp . then by a computation, it can be expressed 



as 



u{t) = S{t)uo = e y ' ^uo{y)dy, (2.1) 

Zm sm(cL'r) Jg2 
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where ;= (—X2,xi) and S{t) is the evolution operator which can be formaUy written 
as S{t) :— e*5('^^*'^^ ) . This formula, which can be deduced from the three-dimensional 
nonlinear Schrodinger equation with a magnetic field discussed in [7, Sec. 9.1] or [1], defines a 
operator S{t), unitary on L^. Note that this formula is valid only for small time, due to the 
singularity formation for the fundamental solution. For this nonlinear Schrodinger equation, 
Strichartz estimates are available. These estimates, mixed time-space estimates, are exactly 
the same as for So{t) — e^*^. Recall the main properties from which such estimates stem. 
The operator S{t) is unitary on L^, ||S'(t)||^2^^2 = 1. In fact, by the Plancherel theorem, we 
have for any (f> ^ L'^ 



wsmw 



OJ 



^ 27r|sinwt| 
27r jsincjtl 

LU 



- cot{ujt)--x -y 



(f>iy)dy 



^\^\^ cat(u,t) 



(t>iy)dy 



\smujt\ 



Lj (xcot(cji) — x^)) 



And for < t ^ the operator is dispersive, with ||'5'(t) H^i^^^o ^ g: N since |sini| ^ ^ \t\ 
for |t| ^ ^. Thus, we can obtain similar Strichartz estimates to the linear Schrodinger operator 
e*l^ by the standard methods (c.f. [11]) provided that only finite time intervals are involved 
(c.f. [6]). 

Proposition 2.1. Let I be an interval contained in [0,tt/2uj]. Then, it holds that 

(1) For any admissible pair (7(23), p) (that is, — 1/2 — 1/p for 2 ^ p < 00), there 
exists Cp such that for any ip € L'^ 

11011^2. (2.2) 

(2) For any admissible pairs (7(^1), Pi) and (7(^2), P2), there exists Cp-^_p2 such that 



[ S{t-s)F{s)ds 



L7(Pl)(/;LPl) 

The above constants are independent of I d [0, 7r/2a;] 
The integral equation reads 



^ ^Pl,P2 II^IIl7(P2)'(/;L''2) " 



(2.3) 



u{t) = S{t)uo — iP I S{t — s) \u\ u{s)As. 
Jo 



(2.4) 



Since the initial data belong to the space S, we naturally need the estimates of VS{t)ip 
and xS{t)(j). In fact, from (|2.ip . we can compute and obtain that 



VS{t)(l) = iujxcoi{ijjt)S{t)(j) ~ iujcot{Lut)S{t){x(f>) + iujS{t){x^ 
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and 



S{t)\7-^(t) = iujx^ coi{u)t)S{t)(t) - iuj cot{ujt)S{t){x^(p) - iujxS{t)(l>, 
S{t)V(j) = iujxcot{LJt)S{t)(f> - iujcot{LJt)S{t){x(j)) + iujx^S{t)(l), 



which yield to 



VS{t)<i)=cos{ijjt)S{t){cos{u)t)V - sm{ujt)W^)(l> 

— iLos\Yi{ijt)S{t) [(cos(a;i)x — sin(ajt)a;"'")(/)] 
xS{t)(j> =cos{uit)S{t) [(cos(u;t)x — sin(a;t)x"'")(/)] 



sm{Ljt)S{t){cos{ujt)V - sm{ujt)\/-^)(j). 



Thus, we have 



S{t){-iW)(p 

= [lo sm{ujt){cos{ujt)x + sin{u;t)x'^) — z cos(a;t)(cos(wi)V + sin(wt)V^)] S{t)(j), 



and 



S{t)wx^ 

= [u cos{u)t){cos{iL)t)x + sm{ujt)x'^) + ism{ut){cos{u)t)W + sin(a;t)V"'")] S{t)(l). 

For convenience of computations, we denote 

J{t) = LO sm{Lot){cos{u)t)x + sin(a;t)a;^) — icos{ujt){cos{Lot)V + sin(a;t)V"'"), 

and the corresponding "orthogonal" operator 

H{t) = to cos{u)t){cos{Lot)x + sin(a;t)a;"'") + i sin(a;t)(cos(a;f)V + sin(a;t)V"'"), 

which will appear in the pseudo-conformal conservation law and play a crucial role 
nonlinear estimates. 
Thus, we get 

J{t) = S{t){-iV)Si-t), H{t) = Sit)LJxSi-t). 
By computation, we can obtain the following commutation relation 



(2.5) 
(2.6) 



(2.7) 

(2.8) 

in the 



1 



. ,2 



J{t),idt + -A- — \xr+LjL, 



= 0, 



1 . w^ ,2 



H{t),tdt + -A- — \x\'+ujL, 



(2.9) 



0. 



In addition, denote M{t) = e"*"^ t'^'iM) and Q{t) = g^'^T-'^otM), then 

J(t) = -icos{Lut)M{t){cos{ujt)V + sin{ujt)V-^)M{-t), 
H{t) =ism{ut)Q{t){cos{u!t)W + sin(wt)V-^)Q(-t). 



(2.10) 
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3 The conserved quantities 

Proposition 3.1. Let u be a solution of the equation p.Sp with the initial data (j> G S(M^). 
Then, we have the following conserved quantities for all t ^ 0: 
(1) The L^-norm: 



II«WII2 = II"0||2 

(2) The energy for the non-rotating part: 



1 



2 , ^ 



(3) The angular momentum expectation: 



/3 



|2ct+2 
l2cr+2 



Eo{uo). 



(3.1) 



(3.2) 



(3.3) 



{L,){t) = / uL,udx= {L,){0). 
Proof. For convenience, we introduce 

eq{u) iut + 2~ ^^^^ u — (3 u + loL^u. 

It is clear that p.ip holds by applying the L^-inner product between eq{u) and u, and 
then taking the imaginary part of the resulting equation. 

Since we can use the identity (|3.3|) in the proof of (|3.2p . we derive p.3p first. Differenti- 
ating (Lz)(t) with respect to t, and integrating by parts, we have 



dt 



/ [ut{x2dx^u - xidx^u) + u{x2dx^ut - xidx^ut)]Ax 



-iut{x2dxiU - xidx^u) - iut{x2dxiU - xidx^u)] Ax 

1 \ 

-Au — \x\^ u- 13 \u\^" u + ajL^uj {x2dxiU - Xidx2u) 



1 

-Au — \xf u- 13 u + ujLzU ) {x2dxiU - xid^^u) 



dx 



Re{Au{x2dxiU - xid^^u)) - -y (Nl^ {x2dxi \uf - xid^^ \u\^)) 
- /3Re(|u|^'^ {x2dxt \u\^ - xid^^ \u\^))dx 



= ^ / {2xiX2 - 2x2Xi \u\^))Ax 

2 JR2 



which yields the desired identity (|3.3p . 

Next, we prove the energy conservation law for the non-rotating part 

Re(eg'(M), ut) = 0, 

where (•, •) denotes the L^-inner product. From the above, we can get 

\dt Wu\' + ^dt \xu\' + -^dt \u\"'+^ + ^dtiuLzu)] dx = 0, 
2 2 (T + i 2 



We consider 



which implies the identity (|3.2p with the help of (|3.3 



□ 
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Remark 3.2. For the equation ()1.5p . the pseudo-conformal type conservation laws are also 
valid: 



2 , 2/3 Sin" (cot) „ 2^+2 , 2/3[aco - 1] 

2 



||g(t)u||- + -:" Ik|l2::2+ ^":: . " / sin2u;s||u(5)||^^+^ci. 



and 



2/3 



\\jit)ur + - \\u\\zi2 = llvz.oii2 + ^ ii«oi 

cr + i cr + i 

2/3[crw-l] /■* 



2tT+2 
20-+2 



3[(JUJ - 1 /■ . „ II , N||2o-+2 , 

4 Nonlinear estimates and the proof of Theorem 11.11 

By computation, we can get, with the help of (|2.10p . that 



J{t) \ur u={a + l) \ur J{t)u - a Wr^' u'j{t)u, 
which implies, in view of^+£=^ + ^ with < e < that 



J{t)\utu ^C\\u\\l^AJ{t)u\\ 



LP 



From the Sobolev embedding theorem and the Holder inequality, it yields 



J(<)|u|'"u 



^(l-pe) (0,T;l( ) ) 

\\J{t)u\\ 



Similarly, we have 



■*( l-pE ) /n T-l( ) "( 



L'Vi-p=7 (0,T;LU-pe; ) 
2cr 



lL°°(0,T;_f/i) 



||i/(t)ul 



Lt(p)(0,T;Lp) ' 



and 



|2cr 
U U 



. ( p V f p V ^ CT 



II"IIl°°(0,T;H1) II"I1l-'(p)(0,T;Lp) 



For convenience, we denote 

III^IIU := \\u\\j, + \\J{t)u\\^ + \\H{t)u\\j,. 
where A denotes a normalized space. Thus, we have 



\u\ U 



■*( i-pe ) Co T-l( i-ps ) 



2cr 



5$CT^ " tCp) ||u||i""oo(o,T;ffi) III"IIIl^<p)(0, T;LP) 



(4.1) 
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For any p £ [2, oo) and M ^ 2C HuoH^, define the workspace {V, d) as 

V := {u : |||u|||^oo(o,T;L2)nL^(p)(0,T;LP) ^ -^}' 

with the distance 

d{u,v) = w||li^(p)(o,T;LP) • 

It is clear that (2?, d) is a Banach space. Let us consider the mapping £^ : {V, d) (V, d) 
defined by 

^ : u{t) ^ S{t)ua -iP [ S{t - s) \u\^" u(s)ds. 
Jo 

For u G {T>,d), by the commutation relation (|2.9p . Proposition 12.11 and the nonlinear 
estimate (|4.ip . we obtain 

^llL°=(0,T;ffl) III^IIIl^(p)(0,T;Lp) 

«CM, (4.2) 

where we have taken T e (0, n/2uj] so smah that Cr^~^^^M2'^ ^ 1/2. Similar to the above, 
a straightforward computation shows that it holds 

d{^ U, ^ V) S^CT^'''~ (l|w|li~(0,T;Hl) + lkllL~(0,T;//l)) III" - «IIIlMp)(0,T;Lp) 

^^d{u,v). (4.3) 

Hence, ^^J' is a contracted mapping from the Banach space (2?, d) to itself. By the Banach 
contraction mapping principle, we know that there exists a unique solution u € (2?, d) to 
p.5p - (|1.6|l . In view of the conservation laws, we can use the standard argument to extend 
it uniquely to a solution at the interval [0, n/2Lu] which satisfies for any t G [0, n/2Lu] and 
P e [2, oo) 

u{t, x), J{t)u{t, x), H{t)u{t, x) G C(0, 7t/2uj; L^iE."^)) n L'^^p^O, Tr/2u; LP{M.'^)). 

Then, we can extend the above solution to a global one by translation. In fact, in order 
to get the solution in the interval (7r/2w, tt/w], we can apply a translation transformation with 
respect to the time variable t such that the initial data u{ti/2lo) are replaced by u{Q). Let 
u{t, x) :~ u{t — 71/210, x), then we have from the original equation with initial data u(7r/2aj, a;) 

1 

lilt + -Au = — |a;|^ u + [3 \u\^" u - ujL^u, a; G M^ t ^ 0, (4.4) 
u(0, x) = uoix) := u(7r/2w, a;), x G (4.5) 

In the same way, we can get a solution u{t,x) of (|4.4p - (|4.5p for t G [0,7r/2w]. It is also a 
solution u{t,x) of (|1.5p - (|1.6p for t G [7r/2a;, tt/o;] and it is unique. Thus, by an induction 
argument with the help of those conserved identities stated in Proposition 13.11 we can obtain 
a global solution u{t,x) of p.5p - p.6p satisfying for any T G (0, oo) 

u{t,x),J{t)u{t,x),H{t)u{t,x) G C(R;L2(]R2)) (0,T;L''(M2)). 

Therefore, we have completed the proof of the main theorem. 
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